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We describe the generation of powerful dispersive waves that are observed when pumping a dual concentric 
core microstructured fiber by means of a sub-nanosecond laser emitting at the wavelength of 1064 nm. The 
presence of three zeros in the dispersion curve, their spectral separation from the pump wavelength, and the 
complex dynamics of solitons originated by the pump pulse break-up, all contribute to boost the amplitude of 
the dispersive wave on the long-wavelength side of the pump. The measured conversion efficiency towards the 
dispersive wave at 1548 nm is as high as 50%. Our experimental analysis of the output spectra is completed by 
the acquisition of the time delays of the different spectral components. Numerical simulations and an analytical 
perturbative analysis identify the central wavelength of the red-shifted pump solitons and the dispersion profile 
of the fiber as the key parameters for determining the efficiency of the dispersive wave generation process. 

OCIS codes: (190.4370) Nonlinear optics, fibers; (190.5530) Pulse propagation and temporal solitons; 

(320.6629) Supercontinuum generation; (060.4005) Microstructured fibers. 
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1. Introduction 

Dispersive waves (DWs), also known as Cherenkov or 
nonsolitonic radiations, are generated in fibers when op¬ 
tical soliton pulses are perturbed by higher-order dis¬ 
persion terms. As initially predicted by Menyuk and 
co-workers for standard fibers, an intense pulse whose in¬ 
put spectrum is close to the zero-dispersion wavelength 
(ZDW) and in the presence of third-order dispersion 
(TOD) sheds power in the form of a long pulse, whose 
spectrum is a narrow peak 0, i- By assuming that 
the TOD and the fourth-order dispersion coefficients are 
known, and by relying on an accurate formula for the 
perturbed soliton, Akhmediev and Karlsson have shown 
that it is possible not only to calculate the DW fre¬ 
quency, but also its temporal shape and the amount of 
radiated energy Q . Even when the full dispersion curve 
is taken into account, the DW wavelength can be read¬ 
ily calculated by imposing a matching condition between 
the phases of the soliton and the DW. In fact it can be 
shown that, depending on the shape of the dispersion 
curve, the DW may be generated at shorter or longer 
wavelengths than that of the pump. Moreover, two or 
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more DWs can be simultaneously generated [ill. 

Microstructured fibers (MFs) allow for a great free¬ 
dom in the desim of the dispersion curve and mode 
effective area P-Q: the former feature is essential to 
control the DW spectral position and the latter helps to 
maximize the nonlinear coefficient. In fact, many groups 
have reported an efficient DW generation (for instance: 
EMI) when the ME is pumped by pulses whose dura¬ 
tion is of a few hundreds femtoseconds or even shorter. 
When sufficiently intense input pulses are used, the DW 
generation is accompanied by the formation of a super¬ 
continuum (SC) spectrum, whose bandwidth can reach 
the extension of a full octave E3“E1| ■ Different phenom¬ 
ena contribute to SC generation: self-phase modulation 
(SPM) and modulation instability (MI) lead to the ini¬ 
tial spectral broadening, which can be further enlarged 
by four-wave mixing (FWM) EMli , soliton fission ( 2 T 
[ 2 ^ and Raman soliton self-frequency shift (2^. More¬ 
over, the large number of pulses fostered by the break¬ 
up of a long pump pulse can also interact with the DW, 
leading to the onset of new spectral peaks and contribut¬ 
ing to determine the extension and flatness of the output 
spectrum [23-1^. 

The development of practical DW sources is gener¬ 
ally limited by the small efficiency of the energy transfer 
from the pump to the DW. In most of experiments re- 
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ported to date, the power carried by the DW is only a 
small percentage of the input pump power; it is possi¬ 
ble to increase the conversion efficiency (up to 65% in 
recent experiments) only by resorting to a femtosecond 
laser source (like a Ti:sapphire laser) [H, [2^ . In 

this case, intense DWs have been observed in the vis¬ 
ible range (even in the violet 0) as well in the near 
and mid-infrared l2i,|3ll. Indeed with ultrashort pump 
pulses their spectrum may be so wide that it incorpo¬ 
rates the resonant condition for DW generation. 

There is a long way ahead in order to achieve high con¬ 
version efficiencies when using nanosecond pulses emit¬ 
ted by low-cost, widely used sources, such as Nd:YAG 
microchip lasers M- Nanosecond pulses cannot directly 
generate DWs, owing to their extremely narrowband 
spectra. However high energy nanosecond pulses may 
break, after a first stage of nonlinear propagation in a 
MF, into a bunch of ultrashort pulses, which in turn can 
generate DWs. 

It must be emphasized that a pulse that sheds light 
into a DW is not necessarily restricted to a fundamental 
or higher-order soliton. In fact, it has been shown that 
intense pulses traveling in the normal dispersion regime 
of the fiber and satisfying a phase-matching relation may 
also effectively build up a DW d, [S^- Moreover, DWs 
have also been observed in a line-defect photonic crys¬ 
tal waveguide with a length of only 1.5 mm: the 
measured 30% conversion efficiency could be explained 
by means of the locking of the velocities of the pump 
soliton and the DW. 

We have recently reported the experimental observa¬ 
tion of a gigantic DW, generated inside a dual concen¬ 
tric core MF pumped by a microchip laser, and carrying 
up to 50% of the input pump power at the fiber out¬ 
put [ 3 ^. In the present work, we further clarify, by 
means of an extensive experimental analysis, the phys¬ 
ical mechanism for such huge energy transfer into the 
DW. Section 2 gives an account of the spectra measured 
when pumping at the wavelengths of 1064 nm and 1030 
nm. Whereas Section 3 describes the spectro-temporal 
analysis of a 200 nm wide bandwidth centered around 
the emitted DW. Section 4 is devoted to the numerical 
and analytical study of the DWs: it also unveils the roles 
of the pump soliton central wavelength and of the fiber 
dispersion profile. Finally, Section 5 briefly summarizes 
the results of our study. 

2. Experiments with two different laser sources 

Our dual concentric core MF has an inner core and a sec¬ 
ond external concentric annular core. These two cores 
are obtained by filling some of the holes in a triangular 
lattice; the holes have a radius of 0.65 fxm and are sepa¬ 
rated by a pitch of 2.6 /im. A scanning electron micro¬ 
scope (SEM) image of the fiber cross section is displayed 
in the inset of Fig. [TJ the two glass hexagonal cores are 
easily recognized. The central core is also doped by Ger¬ 
manium in order to increase both local refractive index 
and nonlinearity [ssl - ls^ . This kind of double core struc¬ 
ture has been proposed and demonstrated to be a very 


effective design to control the magnitude and sign of the 
group velocity dispersion (GVD) in a selected spectral 
range [H, 

The linear guiding properties of our double core MF 
can be calculated from the SEM picture through a nu¬ 
merical mode solver: Fig. [IJa) and Fig. [IJb) show the 
group velocity (l/Pi) and its dispersion (/32), respec¬ 
tively. Looking at the GVD it is interesting to ob¬ 
serve the presence of a large positive dispersion peak: 
as a result, wavelengths in the spectral region around 
1650 nm are expected to overtake spectral components 
at 1450 nm by 31.5 ps/m. This condition is quite un¬ 
usual, since in conventional solid core microstructured 
fibers wavelengths around 1650 nm are in the anomalous 
dispersion region. Hence typical SG generation exhibits 
wavelength components around 1650 nm that appear in 
the trailing edge of the output pulse. According to our 
numerical results, the normal dispersion peak reaches its 
maximum amplitude at 1515 nm, and it is bounded by 
two ZDWs at 1353 nm and at 1669 nm, respectively. It 
is worth noticing that in our fiber dispersion is anoma¬ 
lous in the range between 1018 nm and 1353 nm: as 
a consequence, solitons generated by the break-up of a 
pulse centered at 1064 nm (or at 1030 nm) may only ex¬ 
ist in this wavelength range. The fundamental mode of 
our MF is well confined inside the central core for wave¬ 
lengths shorter than 1400 nm. Whereas for wavelengths 
greater than 1400 nm an important fraction of the op¬ 
tical power is confined by the external hexagonal core 
(as for the fundamental supermode of the dual concen¬ 
tric core fiber discussed in Ref. [1^). It is this abrupt 
change of the mode profile and consequently of its dis¬ 
persion relation which gives rise to the dispersion peak. 
For instance, the calculated effective area at 1064 nm is 
6.2 and it increases up to 41.6 at 1550 nm. 
These numerical results were experimentally confirmed 
by inspecting the mode profile at the fiber output by 
means of an infrared camera. 

By inserting the numerically calculated effective re¬ 
fractive index profile in the well-known soliton-DW 
phase-matching relation (see for instance Refs. Sim, 
[l^[^), the wavelengths of the DW radiation peaks can 
be easily calculated. For our dual core MF, all resonant 
DW wavelengths are plotted as a function of the soliton 
wavelength in Fig. [51 Here the small contribution due 
to the soliton nonlinear phase shift has been neglected. 
We verified the validity of this approximation by evalu¬ 
ating the nonlinear contribution to the DW wavelength: 
for a soliton with a peak power smaller than 10 kW, the 
longer DW wavelength (upper curve of Fig. [5]) increases 
by less than 5 nm. 

In the wavelength range considered in Fig. [5] there 
are always two resonant DWs: for a pump soliton at 
1064 nm, the calculated DW wavelengths are 928 nm 
{DWi, blue curve) and 1587 nm {DW 2 , green curve). 
A noticeable feature is that the spectral position of the 
DWi spans a window which is two times wider than that 
of the mate DW 2 - 
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Fig. 1. Numerically calculated parameters of the double core 
MF: (a) group velocity and (b) its dispersion. The verti¬ 
cal solid red line indicates the 1064 nm wavelength; verti¬ 
cal dashed green lines indicate the three ZDWs (1018 nm, 
1353 nm, 1669 nm). The inset shows a SEM photograph of 
the fiber cross section. 



Soliton wavelength [nm] 

Fig. 2. Calculated positions of DW\ (blue line) and DW 2 
(green line) versus the pump soliton wavelength. The vertical 
lines indicate the wavelengths of 1064 nm (dashed red) and 
1100 nm (solid red). 


In this work, two different laser sources are used as 
the pumps to obtain spectral broadening and DW gen¬ 
eration in the same 4 m long dual core MF. The first 
source is a compact microchip Nd:YAG laser (similar to 
the one used in Ref. ©[II), whose pulses centered at 
1064 nm have a duration of 900 ps. The input beam was 


carefully focused on the fiber facet in order to maximize 
the coupling with the fundamental mode. The second 
source is a mode-locked fiber laser emitting 120 ps pulses 
at 1030 nm: its output fiber pigtail was directly spliced 
to the MF. In both cases, the peak power injected inside 
the MF fiber can be varied up to a maximum value of a 
few kilowatts. 

The spectra obtained by using the microchip laser at 
1064 nm are shown in Fig. [Sj for four different values 
of the peak power injected inside the MF. At the lowest 
power level of 0.1 kW, a broad red-shifted spectral peak 
appears: the maximum of this peak is 17 dB lower than 
the residual pump, and it is located at about 1100 nm; 
this provides a first evidence that the broadening mech¬ 
anism is quite asymmetric around the pump, thus fa¬ 
voring the energy transfer towards longer wavelengths. 
Since this spectrum is measured after 4 m of propaga¬ 
tion, the two MI peaks are only scarcely discernible. 
On the other hand, when measuring the spectra after 
a propagation of less than 1 m the initial growth of two 
almost symmetric MI peaks is apparent (as in the exper¬ 
imental data reported in Ref. |34|). By increasing the 
input pump power, the output spectrum broadens fur¬ 
ther, and a DW (henceforth DW 2 ) peak grows around 
1510 nm. The peak at 1100 nm is still present, and it is 
roughly 5 dB higher than the spectral density dip at the 
long wavelength side of the pump. Note that spectral 
broadening towards shorter wavelengths abruptly stops 
at about 800 nm. At the maximum injected pump power 
of 4 kW, the DW 2 peak is shifted to 1548 nm, with a 
spectral intensity that is only 4 dB lower than the resid¬ 
ual pump peak. Quite strikingly, the DW 2 carries an 
impressive 50% of the total output average power. We 
also underline that the beam image at the fiber output 
(observed with an infrared camera around 1500 nm) did 
agree quite well with the profile predicted by the mode 
solver. 

In the measured spectra, there is always a spectral 
hump around 1100 nm which corresponds to the largest 
fraction of solitons generated by the input pulse break¬ 
up. By considering an average soliton central wavelength 
of 1100 nm, the DW 2 peak is expected to be at 1575 nm 
(see Fig. HD, which is reasonably close to the measured 
value; the residual discrepancy is ascribed to a small 
error in our estimation of the MF refractive index trans¬ 
verse profile. The other resonant dispersive wave peak 
DWi, which is predicted to occur at 874 nm, was not ob¬ 
served in the measured spectra, at least in the form of a 
sharp and isolated peak. In fact the spectrum broadens 
towards shorter wavelengths in a rather uniform way: 
its short-wavelength edge, located at around 800 nm, 
is limited by the condition of group-velocity matching 
with the infrared part of the SC (see Fig. [T](a)). Hence 
its location is only slightly affected by variations in the 
input pump power. We expect that in our experiments 
the DWs are shed by a large number of solitons with a 
distribution of peak powers and time widths: this situa¬ 
tion is markedly different from the case of DWs that are 
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Fig. 3. Measured output spectra when the source is the 
Nd:YAG microchip laser at 1064 nm for four different lev¬ 
els of injected input peak power:(a) 0.1 kW, (b) 0.6 kW, (c) 
2 kW, (d) 4 kW. 


emitted by means of stable femtosecond laser sources. 

Solitons generated by the 1064 nm pump do not 
remain confined around 1100 nm, but travel towards 
longer wavelengths, and give rise to a long spectral tail 
in the 1100-1400 nm range (see Fig. [3tb)): this spec¬ 
tral red-shift is due to TOD, and is enhanced by the 
Raman effect [l^. The more the pulses are red-shifted, 
the more they slow down (see Fig.[TKa)): for this reason, 
solitons interact and frequently overlap in time. Collid¬ 
ing solitons may significantly enhance the efficiency of 
DW generation with respect to the single soliton case, 


see Refs. (40l - l^ . In our experiments, the large number 
of collisions among the solitons originating from pump 
pulse break-up will thus largely contribute to reinforce 
the energy transfer towards the DWs. 

Note that the soliton tunneling mechanism [43l - l45j 
may also lead to a high DW conversion efficiency. How¬ 
ever we may rule out the soliton tunneling mechanism for 
the generation of DW 2 - In fact, the DW 2 spectral peak 
is located well inside the normal dispersion region, and it 
starts growing well before solitons approach the ZDW at 
the border of the normal dispersion barrier (whose calcu¬ 
lated value is 1353 nm). The initial stages of the spectral 
broadening process have been reported in Ref. [^. An 
experimental spectrum similar to the one that we obtain 
at high powers (Fig. [3Kc)-(d)) was previously reported 
for a solid core, double zero dispersion wavelength MF 
by Chapman et al. (4^ . who measured a broad and in¬ 
tense peak at 1.98 /xm beyond a 5 dB flat SC, but in 
that experiment the DW peak only grew after that the 
red-shifted solitons had approached the barrier of nor¬ 
mal dispersion. 

We also observe that our spectra exhibit a hump close 
to 1400 nm (this is clearly visible in Fig. [3Kb)- (c)), which 
could be due to the accumulation of solitons, if we sup¬ 
pose that the actual ZDW in our MF is situated more 
than 50 nm above the numerically estimated value. The 
growth of a DW in the infrared was also observed by 
Kudlinski et al. [l^ by using a MF with two zero dis¬ 
persion wavelengths: however in those experiments the 
DW was about 15 dB lower than the residual pump, 
likely because of the low energy or the limited spectrum 
of the accumulated solitons. 

Once again, it must be underlined that the DW that 
is observed in Fig. |3|is so intense because it is fueled by 
the solitons as soon as they are formed from the break¬ 
up of the pump pulse. Moreover, these solitons carry a 
large fraction of the total energy that is coupled inside 
the MF. The mechanisms that are responsible for the 
exceptional growth of DW 2 at 1.548 as well as for 
the negligible growth of the spectral intensity of DWi 
around 0.9 /xm will be more deeply investigated in the 
following sections. 

For comparison, we report in Fig. [4] the spectra that 
are measured at the MF output by using the fiber laser 
emitting at 1030 nm. As can be seen, even at very 
high input powers a drastic reduction in DW conver¬ 
sion efficiency is observed. At the lowest level of pump 
power that is shown in Fig.jTKa), two small MI peaks on 
both sides of the residual pump can be identified at the 
wavelengths of 1065 nm and 994 nm, respectively. The 
spectral density distribution preserves the same qualita¬ 
tive shape even when the peak pump power is increased 
up to the value of 2.7 kW (Fig. HKb)), but an isolated 
DW peak appears at the wavelength of 1504 nm when 
the injected peak power reaches the value of 5.3 kW 
(Fig.HKc)). At such power the output spectrum extends 
from 828 nm to 1418 nm (when considering the -40 dB 
level). In Fig. HKc) the DW peak amplitude is 39 dB 























5 


lower than the residual pump at 1030 nm, and the DW 
peak rises up to -25 dB at the maximum available power 
of 8 kW (Fig.Hd)). 
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Fig. 4. Measured output spectra when the source is the 
mode-locked fiber laser at 1030 nm for four different levels 
of injected input peak power:(a) 1.37 kW, (b) 2.7 kW, (c) 
5.3 kW, (d) 8 kW. 






ing the fiber laser. We may thus conclude that the DW 
emission efficiency and its spectral position are strongly 
connected with both the carrier wavelengths of the soli- 
tons and the MF dispersion profile. 

3. Spectro-temporal characterization of the DW 

We limit our spectro-temporal analysis to the case where 
a maximum conversion efficiency is obtained, z.e., when 
pumping the MF at 1064 nm. Figure [5] presents the 
experimental setup that has been used to measure the 
relative time-delay among the different spectral compo¬ 
nents of the DW. A cube polarizer (CP) splits the optical 
beam of the microchip laser in two directions in order 
to obtain, by means of a photodiode (A), a trigger sig¬ 
nal to be used in an electrical 16 GHz bandwidth digital 
sampling oscilloscope (DSO). Light is injected into the 
inner core of the 4 m long MF by means of a micro¬ 
lens (L2); the outgoing spectrum is collimated through 
a micro-lens (L3), dispersed by means of a diffractive 
grating (DG) and an aperture of 1 mm: this structure 
works as a tunable bandpass filter with a bandwidth of 
6 nm. Each spectrum slice is then first measured by an 
optical spectrum analyzer (OSA) and then by a photo¬ 
diode (B), which can be interchanged by means of a lin¬ 
ear translation stage. It is thus possible to measure the 
relative time delay between the arrival time of a given 
spectrum slice in B and the trigger signal in A. It is wor¬ 
thy to remind that the microchip laser is operating in 
Q-switching mode, thus the pulse-to-pulse timing jitter 
needs a self-referencing with the time of emission of each 
pulse, and this is the role of photodiode A. Both photo¬ 
diodes A and B are InGaAs PIN diodes, with 12.5 GHz 
bandwidth. 



Fig. 5. Experimental setup for the spectro-temporal charac¬ 
terization of the DWs emitted by the dual concentric core 
MF; L1,L2,L3,L4 are micro-lenses and the A/2 elements are 
half-wavelength plates to control the polarization state (see 
the text for a detailed description of the setup). 


The laser operating at 1030 nm has a pulse duration 
that is about 7 times shorter than that of the microchip 
laser at 1064 nm. Although it is reasonable to suppose 
that a reduction in pulse duration causes a reduction in 
the spectral broadening and hence in the popula¬ 
tion of solitons, with a consequent drop in the emission 
of DWs, this cannot entirely explain the dramatic drop 
of DW generation efficiency which is observed when us- 


We show in Fig. |6Ka) a collection of temporal pro¬ 
files corresponding to different slices of the DW; panels 
(b) and (c) of the same figure report summary graphs 
for the pulse full-width at half-maximum (FWHM) and 
the pulse delay recorded at different central wavelengths. 
The spectral region at 1650 nm leads the wavelengths 
around 1450 nm by 163 ps after 4 m of fiber, which fits 
fairly well with the 31.5 ps/m x 4 m=126 ps of time de- 
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lay predicted by the numerical analysis of the MF guid¬ 
ing properties. However, we underline that the pulse 
duration that is used in our experiment is longer than 
the time delay under test. Therefore our conclusions are 
drawn by estimating the shift of the center of mass of 
the recorded temporal profiles. From our measurements, 
it is confirmed that the infrared part of the spectrum is 
present in the leading edge of the pulse, owing to the 
normal dispersion region of the MF. To conclude this 
section, we may note that the measured time delay re¬ 
mains nearly unchanged for all considered power levels, 
as the delay is related to the linear guiding properties of 
the dual concentric core MF. 
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Fig. 6. Spectro-temporal experimental characterization of 
DW 2 when the pnmp is the Nd:YAG laser: (a) spectro¬ 
gram around the central wavelength of DW 2 ; (b) DW 2 pulse 
FWHM versus wavelength; (c) DW 2 pulse delay versus wave¬ 
length (the choice of the time origin is arbitrary). The filled 
circles in (b) and (c) are the measured values and the con¬ 
tinuous lines are guides for the eye. 


4. Intensity of the DW upon pump wavelength 

In this Section we present a detailed analytical and nu¬ 
merical study of the generation of dispersive waves in 
our dual-core MF, in order to further clarify the phys¬ 
ical mechanisms which may explain the observed large 


conversion efficiencies and their dependence upon the 
pump wavelength and MF dispersion profile. By nu¬ 
merically solving the Generalized Nonlinear Schrddinger 
Equation (GNLSE) by means of the split-step Eourier 
method, the input pulse break-up and the resulting for¬ 
mation of solitons can be calculated both in the time 
and in the frequency domains. In our numerical sim¬ 
ulations we did not include the variation of the mode 
effective area, since it only has a minor contribution on 
the qualitative dynamics of the propagation of solitons 
in the region 1064-1400 nm. 

Let us first consider for simplicity a long square 
shaped pulse with peak power P=1.5 kW and width of 
Ti=500 ps: this ansatz could represent the internal part 
of a sub-nanosecond laser pulse. If we consider a central 
wavelength of 1064 nm, the numerically calculated group 
velocity dispersion is /32 = —5.88 x 10“^^ s^/m. By as¬ 
suming a nonlinear coefficient 7 = 21 x 10 “^ 
we may estimate a MI frequency Jmi = 16.4 THz, 
hence a modulation period Tmi = 61 fs. Erom the ra¬ 
tio Tl/Tmi, we expect to obtain about 8000 solitons, 
with individual energy equal to Es = P/f mi = 91 pJ. 
When considering only first-order solitons, we may cal¬ 
culate a corresponding reference time duration param¬ 
eter To = 2\j32\/{'^Es) = 6 fs, which corresponds to a 
soliton FWHM of 10.6 fs, and peak power Ps = |4lsp = 
£'s/(2To) = 7.7 kW. Although this reasoning is only ap¬ 
proximate (the pulse widths are of the order of just a 
few optical cycles), the previous values may provide a 
first-order estimation of the magnitude and number of 
solitons that are generated from the decay of the initial 
long pump pulse. 

From our full numerical solutions of the GNLSE with 
either single sub-nanosecond input pump pulses or with 
trains of femtosecond solitons (not shown here), we may 
conclude that the formation of the observed huge DW 
peak cannot be explained if the input solitons are in¬ 
dividually considered. In other words, the experimen¬ 
tally observed DW peak cannot be reproduced by sim¬ 
ply adding all the DWs that are generated by each single 
soliton. Actually, when numerically simulating the for¬ 
mation and evolution of a bunch of these solitons along 
the dual concentric core MF, it is apparent that they 
reach different peak amplitudes and central wavelengths. 
Hence, due to their different group velocities, solitons 
interact and undergo multiple collisions. Indeed it was 
already pointed out that soliton-soliton collisions may 
increase the DW intensity by orders of magnitude [iol - 
1^ . Thus we focus our analysis here on the role of the 
soliton parameters in determining the relative amplitude 
of the two predicted DWs: this problem is relevant since 
only one peak (namely, DW 2 ) is clearly visible in our 
measured spectra. 

As discussed in the previous Section, for an input 
pump pulse centered at 1064 nm, during the first stage 
of propagation the generated solitons have wavelengths 
that are mostly located around 1100 nm (Fig.|3Ka)). As 
the measured spectra show, at the end of the 4 m long 
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MF these solitons entirely fill the spectral region where 
dispersion is anomalous (Fig. [31[b)-(d)). Numerical sim¬ 
ulations confirm this picture, where solitons slow down 
and red-shift; the Raman effect gives a contribution to 
this spectral shift, however our numerical results show 
that a considerable red-shift is present even if the Raman 
coefficient is set to zero. 

The blue thick curve of Fig. [7] shows the numerically 
computed spectrum at the output of the 4 m long double 
core MF, computed by assuming an input pump pulse at 
1064 nm, with a duration of 200 ps and a peak power of 
1 kW. The Raman effect was included in our simulations. 
This spectrum can be compared with the measured val¬ 
ues in Fig. ISKd): both the numerical result and the ex¬ 
periment show a residual MI peak around 1100 nm, an 
intense DW 2 and a region of wide spectral broadening 
in the 800-1600 nm range (see also the numerical results 
reported in Ref. 3^). The red thick curve in Fig. [7] is 
obtained by shifting the carrier wavelength of the input 
pump pulse to 1030 nm. The output spectrum is quali¬ 
tatively similar to the previous one, it exhibits a decrease 
in DW 2 , but it does not agree with the experiment of 
Fig. 01 In fact, the pulses emitted by the fiber laser at 
1030 nm are seven times shorter than the pulses emitted 
by the microchip laser at 1064 nm. We thus decreased 
the pump pulse duration to 30 ps: the black thin curve 
of Fig. [7] shows a consistent drop in DW 2 in satisfac¬ 
tory agreement with the experimental data of Fig. 0Kd). 
Additionally, note that the black curve also exhibits a 
sharp sideband peak at 788 nm which is generated by 
four-wave mixing in the early stages of the propagation 
(the associated sideband is at 1486 nm), i.e., before the 
pulse break-up. Fig. HKd) shows the presence of a little 
sideband on the short wavelength limit of the SC, and a 
similar feature is also weakly visible in Fig. [3Kb). 



Fig. 7. (Color online) Numerically calculated output spectra 
(the Raman effect is included) with different input pulses 
with a common peak power of 1 kW. Blue thick curve: pulse 
duration of 200 ps and central wavelength of 1064 nm; red 
thick curve: pulse duration of 200 ps and central wavelength 
of 1030 nm; black thin curve: pulse duration of 30 ps and 
central wavelength of 1030 nm. 

The nanosecond pump pulse may decay in a significant 


number of solitons, nevertheless we focus our attention 
on the propagation of just one of them, since we may 
expect that the soliton central wavelength has an impact 
not only on the DW position, but also on its amplitude. 

As a first example, we show in Fig. [8] the spectral evo¬ 
lution of a 30 fs soliton initially centred at 1100 nm. The 
color code represents the spectral intensity in logarith¬ 
mic scale, and again we included in these simulations the 
Raman effect. Both dispersive waves DWi and DW 2 are 
clearly evident, and the soliton exhibits a marginal red- 
shift. Most of the radiated energy is carried by DWi in 
agreement with the theory that we will illustrate later 
on, however this situation is far from what we observed 
in the experiments. 

Figure |3| shows the evolution of the same 30 fs soli¬ 
ton but with a central wavelength of 1200 nm. Now the 
soliton and DW generation dynamics is completely dif¬ 
ferent from the previous case, and it can be decomposed 
in two steps: a first stage leading to a net red-shift of 
about 30 THz is followed by a relevant emission of DW 2 
that inhibits further red-shift. Surprisingly, such a huge 
soliton self-frequency shift induces a change in the inten¬ 
sity of DW 2 , while leaving its spectral position nearly 
unchanged. 

Larger energy leakages into DWs can substantially re¬ 
duce the soliton red-shift: this is the case of Fig. (TUI 
where the soliton was initially centered at 1300 nm. Fig¬ 
ures ii and [To] show the asymmetry in favour of DW 2 
which is observed in our experiments. More specifically, 
Fig.[n|shows how this effect is intensified after the initial 
stage of red-shift of the soliton central wavelength. 

From the measured spectra of Fig. [3] we may ex¬ 
pect that the pump pulse break-up eventually leads to 
a significant population of red-shifted solitons between 
1100 nm and about 1400 nm. Simulations, such as those 
of Fig. [9] prove that this particular spectral distribution 
of solitons contributes to generate a red-shifted DW 2 
that is much larger than the blue-shifted DWi. 



Fig. 8. Numerically calculated spectral evolution of an input 
soliton of 30 fs centered at 1100 nm. DWi and DW 2 are both 
clearly identifiable. 


It would be interesting to extrapolate the dependence 
of the DW amplitude upon the soliton wavelength, in 
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Fig. 9. Numerically calculated spectral evolution of an input 
soliton of 30 fs centered at 1200 nm. The soliton follows a 
rapid red-shift and then DW 2 is emitted. 



z[m] 


Fig. 10. Numerically calculated spectral evolution of an 
input soliton of 30 fs centered at 1300 nm. There is no no¬ 
ticeable soliton red-shift and the DW 2 is emitted since the 
beginning. 


order to understand the amount of DW 2 that one can 
potentially generate. To this end, we applied the theory 
developed by Akhmediev and Karlsson Q to the specific 
case of our MF fiber. We evaluated the frequencies of 
resonant coupling and the initial DW spectrum for differ¬ 
ent values of the soliton carrier wavelength. To simplify 
the approach we do not take into account the soliton 
self-frequency shifts that originate from either the Ra¬ 
man effect or the high-order dispersion terms. Still we 
do consider the full dispersion profile in order to properly 
localize the dispersive waves. In fact Fig. IH] shows how, 
in a first approximation, red-shift and resonant emission 
can be considered separately. 

Let us start from the GNLSE (with Raman effect now 
neglected) written in the frequency domain as 

BA 

_ _ = 0 ( 1 ) 

where A(t, z) is the complex envelope of the optical field 
and A{uj, z) = -F[A] = A(t, z) exp(iwt) dt is its 


Fourier transform. Equation o is written in a refer¬ 
ence frame moving at the group velocity of the reference 
physical angular frequency ojq. In what follows we set 
the reference frequency (jJq as the carrier angular fre¬ 
quency of a soliton, so that w measures the shift of the 
angular frequency of the optical field uj + uJo from wq. 
The full dispersion profile of the MF is described by the 
function k(uj) = j3{u! + ojo) — /3(oJo) — /5i(wo)w, where 
/3i(wo) has the physical meaning of the reciprocal of the 
group velocity at wq- Note that varying the reference an¬ 
gular frequencies ojq leads to different dispersion profiles 
k(cli). 

The relative time delay between different spectral 
components of the field and the reference frame is pro¬ 
portional to 

f)fzr 

^ ^ =^i(a; + a;o) -/3i(u;o) (2) 

OOJ 

From Eq. @ we can see that k'{uj) measures the group 
velocity mismatch (or difference in group delays) be¬ 
tween the waves at uj + loq and the reference frame at 
Wo- 

In this work, we mainly focus our attention on soli- 
tons that are generated by the long pulse break-up, so 
it is reasonable to represent the linear operator of the 
GNLSE as = /32(wo)w^/2-l-eif(a;), where /32(wo) is 
the group velocity dispersion measured at the reference 
frequency, H{uj) describes higher-order dispersion terms, 
and e is a small parameter that permits to separate short 
and long length scales in the solution of Eq. O- 

We may thus study the solution of Eq. m by using the 
ansatz provided by a first order soliton plus an additional 
small perturbation A(a;, z) = Ao(a;) exp[iKsz]-|-eF(a;, z). 
In our notation, Ao(a;) is the first order soliton solution, 
Ks = 7 |Asp /2 is the nonlinear contribution to the soli¬ 
ton wavenumber, and F{u!,z) is the perturbation. Our 
ansatz neglects the soliton frequency shift that may be 
induced by hence its validity will be limited to the 

early stages of the DW emission process (details on the 
soliton dynamics in presence of third order dispersion 
can be found in Refs. [1, H^-fs^)- Now, since the soliton 
is the first order solution (i.e., with e = 0) of Eq. ([T]), 
we may collect all terms proportional to e and obtain a 
linear forced equation for the perturbation F{ui, z) 

BF 

— - iK(uj)F = iH(uj)A(uj) exp[i/t5z] (3) 

Bz 

The solution of Eq. ©, with zero initial condition, can 
be written as 


F{uj,z) 


H{uj)Ao{lu) 
KS - k{uj) 




(4) 


Equation 0 expresses the well-known property that 
DWs are fed by the spectrum of the soliton, and that 
their amplitude is enhanced at those frequencies ujr 
which satisfy the resonant condition ^(wfl) = ks Q. 
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When dealing with a large number of unequal soil- 
tons, an important issue is the sensitivity of the reso¬ 
nant condition with respect to variations of the soliton 
peak power P 5 . If we perform a series expansion around 
the resonance frequency ujh, we obtain -I- doj^i) ~ 

k{ujr) + dujRK'{ujR) and, by applying a standard error 
analysis for a variation dPs of the soliton peak power, we 
may calculate the sensitivity of the resonant frequency 
to a small change in soliton power with respect to a ref¬ 
erence value 


^ ^ 
dPs 2 k'{u}r) 

Equation ([S]) shows that the variation of the resonant 
frequency resulting from a small variation of the soli¬ 
ton power is determined by the group velocity mismatch 
k'{ujr) = /3i{ujr+uio) — /3i{uio) at the resonant frequency 
ujR + Wo- Large values of the group velocity mismatch, 
which bring long temporal delays between the DW and 
the soliton, will thus result in a small sensitivity of ujr 
to changes of the soliton peak power. 

Equation ([5]) also describes an important property of 
the generated DWs: soliton amplitude jitter may lead 
to fluctuations in the DW frequency, because a change 
in peak power results in a variation of the zeros of the 
phase-matching condition. Note that a similar change 
is also observed if a red-shifting soliton is subject to an 
amplitude reshaping. Since DW frequency fluctuations 
are inversely proportional to k'{ujr), in the presence of 
a distribution of soliton powers, large values of k'{ujr) 
will drastically reduce the standard deviation of the DW 
frequencies. This may explain why the DW 2 peak at 
1548 nm is narrow-band {k!{uir) has a large modulus 
there), while there is virtually no peak but only a flat 
plateau around 800-950 nm for DWi. 

A similar conclusion can also be drawn for the sensi¬ 
tivity of the resonance frequency as a result of changes in 
the soliton central wavelength. From a close inspection 
of Fig. [21 it is clear that the slope of the curve for DWi 
(lower blue curve) around 1064 nm is more than twice 
the slope of the curve for DW 2 (upper green curve). Nu¬ 
merical simulations (not shown here) attest that small 
fluctuations of the soliton wavelength indeed lead to a 
much larger spreading for the DWi resonance than for 
the DW 2 resonance. 

In Fig. |TT]we illustrate the spectrum of the perturba¬ 
tion as it is predicted from Eq. o, for different pump 
soliton wavelengths, and a propagation length of 4 m: 
the generation of DWi, DW 2 and of the red-shifted spec¬ 
tral tail (between 1100 nm and 1400 nm) are clearly rec¬ 
ognizable. Owing to the perturbative nature of F{oj, z), 
the spectrum of Fig. ITT] provides an estimate of the out¬ 
put spectrum at all frequencies, except for values close to 
the soliton reference frequency. This analysis confirms 
that for soliton carrier wavelengths larger than about 
1150 nm the intensity of the DW 2 peak grows larger 
than that associated with DWi. The DW amplitude 
depends not only on the soliton wavelength, but also on 


its FWHM (which is equal to 30 fs in the simulations of 
Fig. llip : this is not surprising since the soliton spectrum 
appears as a multiplying factor in Eq. 0 . 

Since higher-order terms were neglected in deriving 
Eq. 0, we have verified its validity by computing the 
output spectrum through the numerical solution of the 
GNLSE (see Fig. [T^ . The numerically computed DW 
spectra exhibit a good agreement with the theoretically 
predicted spectra of Fig.[TlJ Thus we may conclude that 
the approximate analytical solution of Eq. (jd]) provides 
a quick first estimate of the DW spectrum, which is use¬ 
ful to explain its dependence on the soliton parameters 
and the fiber dispersion profile. Nevertheless, numerical 
simulations are still necessary for an accurate evaluation 
of the DW amplitude, as well as the precise shape of the 
output spectrum. By moving along the different soli¬ 
ton wavelengths of Fig. [Tl]or Fig. [12] we can see how a 
red-shift of the soliton central wavelength can gradually 
intensify DW 2 - 



Fig. 11. Output normalized spectra (in logarithmic scale) 
as given by Eq. 0 versus the frequency detuning a;/27r for 
different central wavelengths of the input first order soliton 
having a FWHM of 30 fs and for 4 m of propagation. The 
region where DW\ [DW 2 ) is generated is highlighted by a 
dashed line (solid line) ellipse. 


5. Conclusions 

We presented an extensive experimental and numerical 
study of the highly efficient generation of a DW spectral 
peak at telecom wavelengths in a dual-core microstruc- 
tured optical fiber, pumped by a near infrared microchip 
laser. Both experiments and simulations agree in their 
predictions that the conversion efficiency into the DW 
spectrum strongly depends on the value of the pump 
wavelength, in combination with the dispersion profile of 
the fiber. Moreover, in spite of the fact that the phase¬ 
matching condition between solitons and DWs provides 
multiple resonances, we found that the conversion effi¬ 
ciency is by far more efficient for the DW that is closest 
to the pump solitons. We further presented an approx¬ 
imate but analytical expression for the DW spectrum 
which explicitly contains the dispersion profile of the 
fiber, and thus provides an useful tool for estimating the 
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Fig. 12. Output normalized spectra (in logarithmic scale) 
versus the frequency detuning a;/27r for different central 
wavelengths of the input first order soliton (FWHM=30 fs) 
calculated by numerically solving the GNLSE for 4 m of prop¬ 
agation. The region where DWi {DW 2 ) is generated is high¬ 
lighted by a dashed line (solid line) ellipse. 


DW position and amplitude. Moreover, the intensity of 
the analytical DW spectrum is proportional to the am¬ 
plitude of the soliton spectral tail at the resonance wave¬ 
length, which confirms the observation that stronger 
DWs are generated when resonance occurs closer to the 
center wavelength of the pumping soliton ensemble. 

We believe that our results can help in understand¬ 
ing the process of highly efficient DW generation in 
dispersion-engineered optical fibers. In particular, our 
analysis may be used as a guideline for the optimization 
and reverse engineering of specialty fibers with the pur¬ 
pose of emitting a target and high energy DW peak in 
any desired spectral region of interest for a particular 
application. 

S.W. is also with Istituto Nazionale di Ottica of the 
Consiglio Nazionale delle Ricerche. We acknowledge 
the partial support from the Region Limousin (C409- 
SPARC), and by the Italian Ministry of University and 
Research (MIUR, Project No. 20I2BFNWZ2) 
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